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Abstract 

In this paper, we prove a Beale-Kato-Majda blow-up criterion in terms of the gradient of 
the velocity only for the strong solution to the 3-D compressible nematic liquid crystal flows 
with nonnegative initial densities. More precisely, the strong solution exists globally if the 
I/ 1 (0, T; L°°)-norm of the gradient of the velocity u is bounded. Our criterion improves the 
recent result of X. Liu and L. Liu f|25|. A blow-up criterion for the compressible liquid crystals 
system. larXiv:1011.4 399v2 [math-ph] 23 Nov. 2010). 
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1 Introduction 

The governing system of equations for the compressible nematic liquid (NLC) crystal flows is 
the following system of scalar or vector fields p(t,x), u(t,x) and d(t,x) for (t, x) G (0, +oo) x f2, 
for a bounded smooth domain C M 3 : 

d t p + div(pu) = 0, 

< dt(pu) +div{pu®u) + VP = [j,Au- AV • (VdOVd- ^(|Vd| 2 + F(d)I)), (1.1) 
d t d+{u-V)d = v{Ad- f{d)) 

together with the initial value conditions: 

p(0, x) = p (x) > 0, u(0, x) = u (x), d(0, x) = d (x), \/x € (1.2) 
and the boundary value conditions: 

u(t,x)=0, d(t,x) =d (x), \d (x)\=l, V(t,x) e [0,+oo) x d£l. (1.3) 
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Here we denote by p, u — (u\, 1*2,1*3), d — (d\, ^2,^3) the unknown density, velocity and orientation 
parameter of liquid crystal, respectively, and P = P{p) is the pressure function. Besides, ji, A and 
v are positive viscosity coefficients. The non-standard term Vd Vii denotes the 3x3 matrix, 
whose (i,j)-th element is given by Y]u— i djdkdjd).- I is the unit matrix. f(d) is a polynomial 
function of d which satisfies f(d) = J^F(cZ), where F(d) is the bulk part of the elastic energy; 
usually we choose F(d) to be the Ginzburg-Landan penalization, i.e., F(d) — jj?(\d\ 2 — l) 2 and 
f(d) — -iz(\d\ 2 — l)d, where a is a positive constant. In what follows, we will assume a = 1 since its 
specific value does not play a special role in our discussion. Besides, we assume that the pressure 
function P satisfies 

P = P(-) e C 1 [0,oo), P(0) = 0. (1.4) 

The above system (jl.ip is a simplified version of Erickscn-Lesile system modeling the flow of 
compressible nematic liquid crystals, and the hydrodynamic theory of liquid crystals was estab- 
lished by Ericksen [6] and Leslie [17] in the 1960's. When d = 0, the system becomes to the 
compressible Navier-Stokes (CNS) equations. Matsumura and Nishida [27J obtained global exis- 
tence of smooth solutions for the initial data is a small perturbation of a non-vacuum equilibrium. 
For the existence of solutions for arbitrary initial value, Lions [18] and Feireisl [9] established the 
global existence of weak solution to the CNS equations. Cho et al.[3][3]|3] proved that the existence 
and uniqueness of local strong solutions of the CNS equations in the case where initial density need 
not to be positive and may vanish in an open set. Xin in |32] showed that there is no global smooth 
solution to the Cauchy problem of the CNS equations with a nontrivial compactly supported initial 
density. Hence, there are many works [3j [TJ |8j HU H3] HH [30] [31] try to establish blow-up criterion 
for the strong solution to the CNS equations. In particular, it is proved in [14] by Huang, Li and 
Xin that the serrin's blow-up criterion (see [55]) for the incompressible Navier-Stokes equations 
still holds for the CNS equations, i.e., if T* is the maximal time of existence strong solution, then 

T !im (II divu|| L i (0i T;L-) + ||p 5 «|U»(o,T ; i,'-)) = 00 (1.5) 

or 

^,S\\P\\^{0,T;L^) + \\P^u\\l^(0.T-L-)) = OO, (1.6) 

where r and s satisfy - + -<l,3<r<oo. In [HHE], Huang et al. established that the Beale- 
Kato-Majda criterion (see [T]) for the ideal incompressible flows still hold for the CNS equations, 
that is 

lim / ||Vu||Loodt = 00. 
t^t*J 

Sun, Wang and Zhang in [30] (see also P3]) obtained another Beale-Kato-Majda criterion in terms 
of the density, i.e., 

lim sup ||/9||l°c(o,T:L°°) = 00. 

T->T* 
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When p is a positive constant, the system (jl.lj) becomes to the incompressible nematic liquid 
crystal (INLC) equations, the global-in-time weak solutions and local-in-time strong solution have 
been studied by Lin and Liu [2UJ IS] • m [H] i Hu and Wang established global existence of strong 
solutions and weak-strong uniqueness for initial data belonging to the Besov spaces of positive 
order under some smallness assumptions. Liu and Cui in 24 obtained that the blow-up criterion 
JO) or (JO) still holds for the solution of the INLC equations. We also refer [TUl [T9j [22j [23j [29] 
and the reference cited therein for other related work on the INLC equations. 

Inspired by the above mentioned works on blow-up criterion of strong solution of CNS and 
INLC equations, particularly the results of Huang etal. [HHE] an d Sun et al. [301131]) we want to 
investigate a similar problem for the compressible nematic liquid crystal flow (|1.1[) - (|1.3[) . Before 
stating the main result, we denote the following simplified notations of Sobolev spaces 

L q := L q (Q), W k ' p := W k ' p (fl), H k := H k (n), := Jfj}(0). 

When the initial vacuum is allowed, the well-posedness and blow-up criterion for strong solu- 
tions to the compressible nematic liquid crystal flows (|l.l|) - (|1.3|l have been investigated by Liu et 
al. in (25] [26] . Here, we write down the main results of Liu et al. [25l [26] . 

Theorem 1.1 Suppose that the initial value (po, uo, do) satisfies the following regularity conditions 

< po e W 1 ' 6 , u G Hi n H 2 and d G H 3 , 
and the compatibility condition 

pAu — A div(Vd Vd — ^(|Vd | 2 + F(do))) - VP(p ) = ^/pg for some g e L 2 . (1.7) 

Then there exist a small T £ (0,oo) and a unique strong solution (p,u,d) to the system with 
initial boundary condition (| 1 - 2[) - (|1.3[) such that 

0<pe C([0, T); W 1 *), p t e C([0, T); L e ), 

ue C{[Q,T); D H 2 ) D L 2 {0,T;W 2 - 6 ), u t G L 2 (0, T; Hq), 
deC{[0,T);H 3 ), d t e C{[Q,T)-Hl)CM 2 {Q : T-H 2 ), 

d tt G L 2 (0, T; L 2 ), y/put G C([0, T); L 2 ). 

Moreover, let T* be the maximal existence time of the solution. If T* < oo, then there holds 

f T 

I™, / (||Vu||£ ( , + ||u|| w ^,«)dt = oo, (1.8) 
where a, f3 satisfying ^ + % < 2 and f3 > 4. 

Remark 1.1 Another similar system of partial differential equations modeling compressible ne- 
matic liquid crystal flows has been studied by Huang, Wang and Wen in jT51[TB] . They obtained the 
existence of local in time strong solution and two blow-up criteria under some suitable assumption 
condition u and d or p and d. 

The purpose of this paper is to obtain the Beale-Kato-Majda blow-up criterion only in terms 
of the gradient of the velocity still holds for the liquid crystal flows. Our main result is the following 
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Theorem 1.2 Assume that (p,u,d) is the strong solution constructed in Theorem ] 1 . 11 and T* be 
the maximal existence time of the solution. If T* < oo, then we have 

limsup||V'u|| L i ( o iT;L oo ) = oo. (1.9) 

T— s-T* 

The proof of this theorem will be given in the next section. As a standard practice, we will 
show that if (jl.9l) does not hold then the strong solution (p, u, d) can be extended beyond the 
time T*. To this end we will step-by-step establish a series of higher-order norm estimates for 
the strong solution (p,u,d). The key fact used in this deduction is that the boundedness of the 
i 1 (0, T; L°°)-norm of Vu implies both boundedness of the L°°(0, T; L°°)-norm of the density p 
and boundedeness of the L°°(0, T; W 1:q )-norm of d with 2 < q < oo. 



2 Proof of Theorem [L2] 

Let (p, it, d) be the unique strong solution to the system (11. ip with initial-boundary condition 
(ll.2l) - (|1.3p . We assume that the opposite to (|1.9[) holds, i.e., 

lim ||V«lUi (0 ,T;i-) < M < oo. (2.1) 

In what follows, we note that C denotes a generic constant depending only on p, A, M, T, Q and 
the initial data. By using the mass conservation equation (ll.ip i and the assumption (I2.ip . it is 
easy to obtain the L°°-norm bounds of the density, 

Lemma 2.1 Assume that 

I ||divtt||ioodt< C, 0<T<T*, (2.2) 

then 

IMU~(o,t ; l«) < C VO<T<T*. (2.3) 

Proof. The proof is essentially due to Huang and Xin [T2], for reader's convenience, we sketch it 
here. 

Multiplying the mass conservation equation (|l.ip i by qp q ~ Y with q > 1, it follows that 

d t (p q ) + div(p q u) + (q- l)p q divu = 0. 
Integrating the above equality over SI yields 

fitllPllI, <(g-l)||div U |M| P ||I 9 , 

i.e., 

atiHiw^^^Hdivuiu-iipiu.. (2.4) 

q 
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The condition (|2.2j) and the estimate (|2.3p imply that 

8t\\p\\L*<C forVg>l, 

where C is a positive constant independent of q, letting q — > oo, we obtain (|2.3I) . and this completes 
the proof of the lemma. □ 

According to the assumption (|1.4p on the pressure P and Lemma |2~TI it is easy to obtain 

sup {\\P(p)\\l~,\\P'(p)\\l°°} < C < oo. (2.5) 

0<t<T 

Now, let us derive the stand energy inequality. 
Lemma 2.2 There holds 



sup / (p\u\ 2 + |Vd| 2 + 2F(d))dx + [ [ \\7u\ 2 dxdt + [ [ \Ad - 
o<t<TJn Jo Jn Jo Jn 



f{d)\ 2 dxdt <C. (2.6) 



Proof. Multiplying the momentum equation (ll.ll) 2 by m, integrating over Q and making use of 
the mass conversation equation (| 1 . 1 [) 1 , it follows that 

\~T f P\ u \ 2dx + P- [ |Vu| 2 dx = - / uVPdx~\ [ (u-V)d-(Ad-f(d))dx, (2.7) 
2 dt J n J u J n J n 

where we have used the fact that div(Vd Vd) = (\7d) T Ad - V ■ Multiplying the liquid 

crystal equation (|l.ip 3 by Ad — f(d) and integrating over O, we obtain 

— / {hvd\ 2 +F(d))dx + v [ \Ad- f(d)\ 2 dx= [ (u ■ V)d • (Ad - f(d))dx. (2.8) 
dt Jn 2 Jq Jn 

Combining (j2~7)) and ([2~8)l together 

— / [i(p|u| 2 + A|Vd| 2 ) + AF(d)]d2; + M / |Vu| 2 da; + \v [ \Ad-f(d)\ 2 dx 
dt Jn 2 Jn Jn 

= - uVPdx= \ Pdivudx<e / iVupdx + Ce^ 1 , (2.9) 
Jn Jn Jn 

where we have used the estimates (|2.3[) . (|2.5p and the Young inequality. Taking e small enough 
and applying the Gronwall's inequality, we can establish the estimate (|2.6p immediately. □ 

In the next lemma, we will derive some estimates of d. 
Lemma 2.3 Under the assumption (|2.ip . it holds that for < T < T* 



sup (\\d\\ Lq + ||Vd|| L9 ) < C for all 2 < q < oo; (2.10) 

0<t<T 

sup ||Vd||ia+ / \\d t \\ 2 L2 dt <C; (2.11) 

0<t<T Jo 



Proof. We first multiplying the liquid crystal equation (|l.ip 3 by q\d\ q 2 d with q > 2, and inte- 
grating over f2, then there holds 



{qv\Vd\ 2 \d\ 2 + q(q - 2)is\d\ q ~ 2 | V\d\ \ 2 )dx 

2 

= -V I Uidi(\d\ g )dx - qv f \d\ q+2 dx + qv f \d\ q dx 
■ , Jn Jn Jn 



i=l 



= -V / diUi\d\ q dx - qv I \d\ q+2 dx + qv [ \d\ q dx 
i=1 Jn Jn Jn 

<a(||vt»iu- + i)||d||£,. 

By using the Gronwall's inequality, one obtains the inequality 

sup \\d\\ Lq < C for all q > 2. (2.12) 

0<t<T 

By letting q — > oo, we notice that the estimate ([2.12j) still holds. 

Multiplying the gradient of the liquid crystal equation (jl.ip , by q\ Vd\ q ~ 2 Vd with q > 2, and 
integrating over f2, then there holds 

d 



dt 

3 



l|Vd||I,+ / MV(Vd)| 2 |Vdr 2 +q{q-2)v\V\Vd\\ 2 \Vd\ q - 2 )dx 
Jn 

3 3 

V f u l( 9 4 (|Vd| 9 )dx- Vg /" Viiidid|Vd|«- 2 Vddz-z/g / V(|d| 2 d)|Vd|«~ 2 Vddx+w? / \\7d\ q dx 

3 3 

=5~] / <9;U;|Vd| 9 d:r-Vg / Vu i a i d|Vci|^ 2 Vddx-^ / V(|d| 2 d)|Vd|«- 2 Vddx+^ f \Vd\ q dx 
i=1 Jn i=1 Jn Jn Jn 

<C||V«||£»||Vd||J, + i>?||Vd||£, —vq [ V(|d| 2 d)|Vd| 9 ~ 2 Vdda; 

Jn 

=(C||V«||£ao+i/g)||Vd||J,,-i/g I \d\ 2 Vd\Vd\ q - 2 Vddx-vq [ dV(|d| 2 )|Vd|^ 2 Vddx 

Jn Jn 

=(C||Vu|| £ ~+i/g)||Vd||£,-3i/g / |d| 2 |Vd| 9 da; 



<C(||Vu|| L »+l)||Vd||i„ 

where we have used the fact that V|d| 2 = 2|d|V|d| = 2|d|^j^ = 2dVd in the last equality. By 
using the Gronwall's inequality again, we obtain 

sup ||Vd|| L , < C for all q > 2. (2.13) 

0<t<T 

Letting q — > oo, estimate (|2.13j) still holds, and the inequalities f)2 . 12[) and (|2.13|) imply that 
estimate (|2 . 10|) holds. 

To prove the estimate (|2.11[) . we multiplying the liquid crystal equation (| 1 . 1 j) 3 by d t and inte- 
grating over Q, then 

||dt||l a + ~ / |Vd| 2 dx = - / (u-V)ddtdx-v [ f(d)d t dx 
z at Jn Jn Jn 
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<c(NU 2 ||w|MWU= + \\4U\\d\\ L 2\\d t \\ L 2 + \\d\\ L 4d t \\ L i) 



<^IKIl! 2 + c, 



where we have used the estimates (|2.6p and (|2.10[) . Integrating the above inequality over [0,T] 
gives the estimate (|2.1ip . □ 

For function /eflx(0,T), let 

/ = /t + u ■ V/ 

denote the material derivative of the function /. Then we have following lemma. 



Lemma 2.4 Under the assumption (|2.1[) , it holds that for < T < T* 



sup (llVullla+HVpllia + \\d\\ 2 H2 ) + / / (p\u\ 2 + \\7d t \ 2 )dxdt < C; (2.14) 
o<t<T Jo Jo 

[ HVdlllra* < <7. (2.15) 
Jo 

Proof. Noticing that the momentum equation (|l.ip 2 can be rewrote as 

pit + VP = pAu - \(Vd) T (Ad - /(d)). (2.16) 

Multiplying the equation (|2.16[) above by u and integrating over 17, one obtains the equality 
p d 



2 dt 



l|Vu||x,2 + / p\u\ 2 dx 



n 

=fi / u • VitAudx + / PdivM t dx— / u ■ VuVPda; 
Jo Jo Jo 

-A [ (ut- V)d(Ad - f(d))dx -A f (u- V)u ■ Vd(Ad - /(d)) da; (2.17) 
Jo Jo 

Combining the mass conservation equation (|l.ip i and the assumption (|2.ip , it follows that the 
pressure P satisfies the following equation 

P t + P'(p)Vp • u + P'(p)pdivu = 0. (2.18) 

Hence, we have 

Pdivu 4 da;=— f P div udx — [ P t divudx 

_ d 
~~dt 



/o Jo 
To estimate the term —A f n (ut ■ V)d(Ad — /(d))dx, we have 



/ P div udx + / P'(p)(Vp • u + pdivu) div udx. (2-19) 
Jo Jo 



-A / (ut ■ V)d(Ad — f(d))dx = A 7^ ( / dju it diddjddx + / uudidjddjddx) + A / u t -Vdf{d)dx 
Jo , in Jo Jo 



=A^(/ djuudiddjddx-^ f d iUu \djd\ 2 dx) - [ dmti^- - ^-)dx 

i j=X i— 1 ^ ~ 

=A^^{— / (djUididdjd — —diUi\djd\ 2 )dx — j djUidid t djddx — j djUididdjdtdx 

+ J diUidjdtfyddx-j j { ^Ml_^pl )dx+ [ dlUl (\d\*d t -\d\d t )dx} 

d F 1 1 1 

<A — / {djUididdjd - -diUi\djd\ 2 - -diU^df + -diUi\d\ 2 )dx 
at Jn £ 4 2 



+ C||Vu|| L2 ||Vd|| L =c||Vd 4 || L 2 + C7||Vu|U=||rf*IU=(||rf||x,«» + 



d f 1 1 1 

<A ^ — / (djUidtddjd- -diU^djd] 2 - -<9 4 w 4 |d| 4 + -5 l u 4 |d| 2 )dx 

+ Ge^HVullia + e||Vdf||i a + C||Vtt||^ + C||d t ||| 2 , (2.20) 



where we have used estimate (|2.10l) in the last inequality. Inserting (|2 . 19[) and (|2.20[) into (I2.17p . 
and integrating over [0, T] give that 

l|Vu||i»+ I I p\ii\ 2 dxdt 
Jo Jn 

<C+c[ [u- \7uAudxdt + C ( P(p)dwudx{T) + C [ [ P'(p){Vp ■ u + pdivu) divudxdt 
Jo Jn Ja Jo Jo. 

3 

+ C J2[ (djUididdjd -]-diUi\djd\ 2 - ^d lUl \d\ 4 + ~diUi\d\ 2 )dx{T) 
ij=i Jn - 

+ e f \\Vd t f L idt + Cs- 1 [ \\Vu\\ 2 L2 dt + C [ (||Vu|||a + ||dt||i>)d« 
Jo Jo Jo 

+ [ f \u\\S7u\\\7P\dxdt + C f f |u||Vu||Vd|(|Ad| + |/(d)|)cLcdt 
Jo Jn Jo Jn 

<C + e [ ||Vcf t ||| 2 dt + C [ [ u-\7uAudxdt + C [ P(p)dWudx(T) 
Jo Jo Jn Jn 

+ C [ f P'(p)(Vp-u + pdivu)divudxdt 
Jo Jn 

3 

+ C J2[ (djUididdjd- ~d iUi \djd\ 2 - \d lUl \d\ A + ^d lUl \d\ 2 )dx(T) 

i,j=l Jn 

+ [ [ \u\\Vu\\VP\dxdt + C [ [ \u\\Vu\\Vd\(\Ad\ + \f(d)\)dxdt, (2.21) 
Jo Jn Jo Jn 

where we have used the estimate (|2.6|) and p. lip . To estimate the terms on the right side of (|2.21|) . 
by using Lemma |2"7T1 the estimates (12.61) . (|2.10p and (|2.11l) . we get 

nu ■ SIuAudxdt = / / (—djUidiudjU — Uididjudju)dxdt 
a i a— i Jo Jn 



= / / (— d-jUidiudjU + -diUi\dju\ 2 )dxdt 
iJ=1 Jo J a 2 

<C [ \\Vu\\ 3 L3 dt < C [ ||Vit||i«»||Vu||| a d* (2.22) 
Jo Jo 

f P(p)dWudx(T)< h\Vu\\\,+C [ \P{ P )\ 2 dx<-\\Vu\\ 2 L2 +C- (2.23) 
Jn 4 Jn 4 

/ / P'(p){\7p-u)dwudxdt < C I HVpll^llull^lldiviillioadt 
Jo Jn Jo 

<C [ T \\Vp\\ 2 L2 \\Vu\\ L ~dt + C [ T \\u\\ 2 L2 \\Vu\\ L ~dt 
Jo Jo 

<C I ||Vp||i a ||Vu||io.dt + C f ||Vu|| L oodt 
Jo Jo 

<C I ||Vp||| 2 ||Vu||ioodt + C; (2.24) 
Jo 

[ [ P'(p)p\divu\ 2 dxd<C [ \\p\\ L oo || Vu|| L2 d< < C [ ||Vu||| 2 di < C; (2.25) 
Jo Jn Jo Jo 

3 r 111 
^2 J (djUididdjd - -diUi\djd\ 2 - -diU^d^ + -d i u l \d\ 2 )dx(T) 



< G(||Vu|| La ||Vd|| £a ||Vd|| L - + ||Vu|U a ||d|| La (||d||i« + ||d|U»)) 

< Jllv u || 2 L2 + cnwi&iivdiii, + c\\d\\U\\d\\ 6 LX + HdHi.) 

<\\\Vu\\ 2 L2 + C (2.26) 



( [ \u\\Vu\\VP\dxdt<C ( [ \u\\Wu\\Wp\dxdt 
Jo Jn Jo Jn 



<C \\u\\ L 4Vu\\ L3 \\V P \\ L ,dt<C ||Vtt||£ 3 ||Vu||2 08 ||Vp||L»d* 
Jo Jo 

< C [ T (\\Vu\\ 2 L2 \\Vp\\ 2 L2 + ||V«||ioo||V«||| a )di 



o 



<C / T (||V U ||i 2 ||Vp||i 2 + ||V U || L »||V U ||i 2 +C7)di; (2.27) 
Jo 

[ f \u\\Wu\\Wd\(\Ad\ + \f(d)\)dxdt<C f (||«|Ua||V«|U6||Vd|Ue||Ad||^ 
Jo Jn Jo 

+ Nl^llVull^HVdll^lldlUco + ||d||£«,))dt 

<C [ (||V M || L2 ||V 2 u|| i2 ||Ad|| L2 + ||V M ||i 2 )di 
Jo 

< [ (e||V 2 ii||| 2 +C£" 1 ||Vw||i 2 ||Ad|| 2 :2 )dt + C. (2.28) 
Jo 

By using the stand elliptic regularity result to (|2.16[) . we have 

||V 2 U || 2 2 <||A W || 2 2 + ||V U || 2 2 < C(||Vu|| 2 2 + \\pu\\ 2 L2 + ||VP|| 2 2 + ||(Vdf(Ad- /(d))|| 2 2 ) 



<c(||v«||£ a + + IIVplH, + ||Vd||i„(||Ad||i a + \\M\\h)) 

<C(\\Wu\\l 2 + \\piu\\l 2 + ||V P ||| 2 + ||Ad|| 2 a + 1). (2.29) 
Combining estimates (|2.21l) - (|2.29p and taking e small enough, we can get 

l|Vu||£a+ / / p\ii\ 2 dxdt 
Jo Jn 

<C + ef \\\7d t \\ 2 L2 dt + C f (||Vp||| a + ||Vu||| a +||Ad||£ a )(l|Vu||£«. + ||Vu||£ a + l)dt (2.30) 



To estimate the orientation parameter d, by the standard elliptic regularity result to the liquid 
crystal equation (II .![) ■-, . one obtains that 

||V 3 d|| L 2 <C(\\Vdt\\ L Z + ||V(« • Vd)|| L 2 + \\Vf(d)\\ L 2 + \\d \\ H3 ) 

<CQ\Vdt\\v + ||V«|| La ||Vd|| L =. + ||| W ||V 2 d||| i2 + ||Vd|| La (||d||io.+||d||L-) + ||d || ff 3) 
<C*(||Vd 4 || L 2 + \\Vu\\ L * + |||w||V 2 d||| ia + C) (2.31) 

Multiplying the liquid crystal equation (ll.l[) o by Ad t , and integrating over f2, then we have 

d 
dt 



f ^|Ad| 2 d2;+ f |Vd t | 2 dx 
Jn Jn 

/ u-VdAd t dx + L< / (|d| 2 - l)dAd t dx 
Jn Jn 



V f Uididdjdtdx-v [ V(|d| 2 d)Vd t dx + v [ VdVd t dx 
■1 Jn Jn Jn 

yj / djUididdjd t dx — / Uididjddjd t dx — v I V(|d| 2 d)Vd t dir + v / VdVd t dx 
n Jn , Jn Jn Jn 



<C(\\Vu\\ L 4Vd\\ L <*>\\Vd t \\ L2 

+ ||wV 2 d|| i 2||Vd 4 || L 2+ ||d|| 2 00 ||Vd|| L 2||Vd t || i 2+ ||Vd|| L 2||Vd t || i 2) 

<e\\Vd t \\ 2 L2 +Ce- 1 (\\Vu\\ 2 L2 + [ |u| 2 |V 2 d| 2 dx + 1), (2.32) 

Jn 

where we have used the Holder inequality and estimates (|2.10l) . For the term J n |m| 2 | V 2 d| 2 da;, 
applying estimate (I2.3ip . we have for 77 > 

f M 2 |V 2 d| 2 da; < C||tt|||,B||V 2 d||ia < C||VM|||2||Vd|| L 6||V 3 d|| L 2 
<r ? ||V 3 d||| a +Cr 1 ||V U ||| a 

<V\W<h\\l^+V [ luflV'dfdx + vWVuWl.+Cv-'iWVuWi.+C). 
Jn 

Hence, taking r\ small enough 

M 2 |V 2 d| 2 dx <2r,||Vd i || 2 2 + C^dlV^lH.dlVwII 2 . + 1) + C). (2.33) 
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Inserting (12.33)) into (|2.32[) . taking e, r\ small enough and integrating above inequality over (0; T] 
ensure that 

[ [ \Vd t \ 2 dxdt<C [ ||V«||ia(||V«||ia + l)dt + C. (2.34) 
Jo Jn Jo 

Now, we will estimate the density p. Applying the operator V to the mass conservation equation 
(11.11) -p then multiplying it by Vp and integrating over fl yield 



d_ 

df 



\\Vp\\ 2 L 2=- / |Vp| 2 divudx - 2 / pVpV divudx - 2 / (Vp • Vu)Vpdz 
Jn Jn Jn 

<C||Vp||i»||Vu||i- +C||Vp|| i2 ||Vdivu|| i2 

<e||V 2 U ||| 2 + Ce^WVpWliiWVuho* + 1) 

<e(\\piu\\h + ||V«||i a + ||Vp||| 2 + ||Arf|| 2 L2 + 1) + ae- 1 ||Vp||i a (||V«)U- + 1) 
<e\\piu\\h + Ce^WVpWUWVuU- + 1) + C(\\Vu\\l* + ||Ad||£ 2 ), 

where we have used the estimate (|2.29|) in the above inequality. Integrating the above estimate 
over (0, T] gives that 



llVpHi, < e / \\piu\\Udt + / (Ce-^lVpHiidlVtilU- + 1) + C(\\Vu\\i? + ||Ad||£ 2 ))dt (2.35) 
Jo Jo 

Combining estimates (|2.30l) . (I2.34|) and (|2.35[) . and taking e small enough, one obtains that 

\\Vu\\l 2 + ||Vp||| 2 + \\Ad\\l- 2 + [ T [ (p\u\ 2 + |Vd t | 2 )dzdt 

Jo Jn 

<C + C [ (||Vp||| 2 + ||Vu||| a + ||Ad|| 2 2 )(||Vw|| L =o + \\Vuf L2 + l)d*. (2.36) 



Since the energy estimate (|2.6[) implies that J Q || Vu|| 2 2 dt < C. By using the Gronwall's inequality, 
the elliptic regularity result ||V 2 <i||i2 < C(||Ad||^2 + ||<io||jj 2 ) and noticing that the assumption 
(|2.1j) . we deduce that the inequality (12. 14)) holds. 

To prove the estimate (12.15)) . by using the standard elliptic regularity result on (jl.lL . we have 

||V 2 d|| 2 3 <C(||d t || 2 3 + \\U ■ Vd|| 2 3 + ||/(d)|| 2 3 + \\d \\ 2 H3 ) 

<C(\\d t \\ L 4Vd t \\ L , + Hli.||Vd||i. + lldllisdldlll- + ||d||| a ) + C) 

<C(\\d t \\l 2 + \\Vd t \\\ 2 + \\Vu\\l 2+ C) 

<C(\\d t \\ 2 L2 + \\Vd t \\l 2 +C), (2.37) 

where we have used the estimate (|2.14l) in the last inequality. Then by using the estimates (|2.10[) . 
(|2.1ip . (|2.14[) and the above inequality, we have 



[ \\Vdf H2 dt< / T (||V 3 d||| 2 + ||Vd||i 2 )dt 
Jo Jo 

<C A||Vd f || 2 2 + ||V(u • Vd)||£ a + ||V/(d)||l 2 + C)dt 
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<c 



' f (\\vd t \\l 2 + ||v«||i a ||vd||i- + ||«||i.||v 2 d||i3 + ||Vd||i a (||d||i« + ||d||i„) + c)dt 

Jo 

<c [ T (\\vd t \\l 2 + \\v*d\\l 2 + c)dt 

Jo 

<C [ (HVdtllii + Hdtllia + C)dt < C. 
Jo 

This completes the proof of Lemma 12.41 □ 
Lemma 2.5 Under the assumption (|2.1j) . it holds that for < T < T* 

sup (||p^ t |l| 2 + ||Vd t || 2 L2 )+ [ T (\\Vu t \\l 2 + \\(Ad-f(d)) t \\l 2 )dt<C. (2.38) 

Proof. Differentiating the momentum equation (|l.l[) n with respect to time, multiplying the re- 
sulting equation by u tl integrating it over f2 and making use of the mass conservation equation 
(jf.ip 1 , one obtains that 



ld_ 

Ydt 



I p\u t \ 2 dx + fi / |Vut| 2 dx— / P t divutdx 
Jo Ja Jn 



n 

|2 



■hi 



- pu- V( 1 ^- + (it • V)uu t ) + p(u t ■ V)uu t dx - X (u t - V)d t (Ad - f(d))dt 



n 2 ./,, 

A / (ut-V)d(Ad- f(d)) t dx 
Jn 

12 l„,.|2 



= / V/9 • u-L^ — |- pdivuJ— ^ pu ■ V((n • V)uut) — p{ut ■ V)uutdx 

Jn 2 2 

-X [ {u t - V)d t (Ad - /(d))d< -X [ {ut- V)d(Ad - /(d)) t dz. (2.39) 
Differentiating the liquid crystal equation (jf ,f L with respect to time gives 
(ut • V)d = i/(Ad - /(d)) t - d tt - (u • V)d t . 
Multiplying the above equality with (Ad — f(d)) t and integrating over f2, one obtains the equality 

(ut-V)d(Ad- f(d)) t dx 

n 

H(Ad - f(d)) t \ 2 - d tt Adt + d tt f(d) t - (u • V)d t (Ad - f(d)) t )dx 



n 



= J v\(Ad - }{d)) t \ 2 dx + ^j t \\Wdt\\h - j (("f ' V)d)/(d) t dz 
+ / vf{d) t {Ad- f{d)) t dx- j {{u ■ V)d t )Ad t dx 



= f v\{Ad- f{d)) t \ 2 dx + ~\\Vd t \\ 2 L2 - f {{u t ■ V)d)f(d) t dx 

+ j uf(d) t (Ad- f(d)) t dx+ [ ((Vu ■ V)d t Vd t - - divu|Vd t | 2 )dx, (2.40) 
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where we have used the fact 



— / ((« • V)dt)Adfdx = — 2, / Uidid t djjdtdx 
Jn iJ=1 Jn 

3 „ 

= z2 {djUididtdjdt+Uidii "^"' )dx 
„• ZTi Jo, 



\dA\ 



3 



^ ( / (djUididtdjdtdx - - / diUi\d 3 d t \ 2 dx) 
I ((Vu-V)dtVd t - idivu|Vd f | 2 )dx 



in the last equality. 

From the equation ()2 . 1 8|) . we can derive 



/ P t div u t dx = - / P'(p)(Vp ■ u + pdivu)dwu t dx. (2-41) 
Jn Jn 

Inserting the equalities (|2.40l) and f|2.41[) into (|2.39j) derives 



dt J n { 2 plutl + A|Wtl )cb + /i||Vu * 11 ^ 2 + v|l(Ad ~ K d Mh 

2 i „l J:, r „,ll„, |2 i „|„,||„, ||Y7„,|2 , „|„,|2|„, ||Y72„,| , „|„,|2| 



<C/ (|V/»|]«]|i*tr+pl div«||«tr+/>|ix||ttt|]V«|^+p|«r|t*t]|V J «|+p|«rlV«||Vi^|+p|titr|V«|)da: 
+ C f(\(u t - V)df(d) t \ + \(Ad- f(d)) t f(d) t \ + |(V« • V)d f Vd d | + | divu||Vd*| 2 )dx 
+ C / |(u t .V)d t (Ad-/(d))|da: + C7 / |p / (p)||Vp||«||divu t |+p|i y (p)||divu||div« t |da: 



3=1 

We will estimate Jj term by term. In the following calculations, we will make extensive use of 
Sobolev embedding, Holder inequality, Lemmas I2.1M2. 41 and the estimate (|2.5I) . 

Ji < C||Vp|| L »||«t||i.||u|U« < C\\Vu t \\ 2 L2 \\Vu\\ L 2 < e\\Vu t \\ 2 L2 + Ce~ x \ 

J 2 + J e < C\\Vu\\ L c\\p?u t \\ 2 L2 ; 

Ji < C7||i**|U«||Vd||xa(||d||ioo + l)||dt|Ua 

< C\\piu t \\ L 2\\Vd t \\ L 2 < e\\piu t \\ 2 L2 + Ce-'WVdtWl.; 

J 8 < C\\(Ad - f(dj) t \\ L2 \\f(d) t \\ L 2 < s\\(Ad-f(d)) t \\ 2 L2 + Ce^WdtW 

,h + Jw= / |(Vli • V)d t Vd t \ + |divu||Vd t | 2 dx 
Jn 

<C||V U || L ~||Vd t || 2 2 ; 
Ju < C , ||u t || i e||Vdt|U»||Ad|| ii 3 + C\\u t \\ L 2\\Vd t \\ L 2(\\d\\ 2 LX + l)||d||io 

< C\\Vu t \\ L2 \\Vd t \\ L2 \\d\\l 2 \\Vd\\l 2 + C\\ P K t \\ L2 \\Vd t \\ L2 



2 . 

L- - 
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< e ||V«t||i a + Ce^WVdtWUmW 2 ^ + 1) + e\\p*u t \\l. + Ce^\\ Vd t ||£ 2 

< e||V«t||i a +e\\p^u t \\l, + Ce-'WVdtWUWVdWl, + 1); 

and 

Jia < C\\p\\ L ~\\Vu\\ L 2\\Vu t \\ L 2 < e\\Vu t \\ 2 L2 + Ce- 1 . 

To estimate the terms J3, J4, J5 and J12, by using the standard elliptic estimate on f|2 . 16[) and 
making use of the liquid crystal equation yield that 

His* <C(\\pu t \\ L 2 + \\pu- Vu\\ L 2 + ||VP|U» + ||(Vd) T (Ad - /(d))IU») 

<C(|| / 9M i || L 2 + \\p U ■ Vm|| L 2 + ||VF|| L 2 + \\(Wd) T {d t + (u ■ V)d)|| L2 ) 
<C(\\p*ut\\w + ||p|M|u|M|Vu|| L « + l|Vp|| L 2 + \\(Vdf(d t + (u ■ V)d)\\ L 2 
<C{\\p^u t \\ L 2 + <t\\Vu\\ h i + a^WVuW^ + \\Vp\\ L 2 + ||Vd|| i3 ||d t || L 6 + ||Vrf||| 2 ||u||i6) 
<C{\\p^u t \\ L 2 + ct||V«||hi + (t _1 ||Vu|U 2 + \\Vp\\v + \\Vdt\\v), 

where we have used the estimates (|2.6p and (|2.10p in the last inequality. Taking <t small enough 
yields 

\\u\\ H * <C(\\piu t \\ L 2 + \\Vu\\ L 2 + \\Vp\\ L 2 + ||Vdt||£»). (2.43) 
Making use of estimates (|2.14[) and (|2.43p . we can estimate J3, J4, J5 and J12 as 

J 3 + J 4 + J5= / p\u\\Vu\ 2 \u t \+ p\u\ 2 \V 2 u\\u t \+ p\u\ 2 \Vu\\Vu t \dx 
Jn 

< C(||p||L-.||«||L«||«t|U«||V«||£a||Vu|| L a + L ~ || u t \\ L e \\ u\\ \ e || V 2 u|| L 2 

+ IIHUHMl!<3||V'it||L6||Vitt|| L 2) 

< C||V«t|U» II V«||| 3 Htill^ < e\\Vu t \\ 2 L2 + Ce- x \\uf H2 

< e\\Vu t \\ 2 L2 + ae- 1 (||p*« t ||i a + ||Vd t || 2 L2 + ||V U || 2 L2 + ||Vp||| 2 ) 

< e||V Ut ||| 2 + Ce-^WpiutWh + ||Vd t || 2 2 + 1); 
J12 < C||Vp||L>||u||i-||V«t||£> < C\\Vu t \\ L 2\\u\\ H 2 

<£\Wu t \\ 2 L2 +Ce- 1 \\u\\ 2 H2 

< e\\Vu t \\ 2 L 2 + Ce-^utWh + \\Vdt\\v> + + l|Vp||£0 

< ellVtttH^ + Ce-Wlp^utWl. + HVdtHia + 1) 

Substituting all the estimates of Jj into (|2.42l) . and taking e small enough, we obtain 

j t J^p\u t \ 2 + IVdtftds + ||V U4 || 2 2 + \\(Ad-f(d)) t \\ 2 L 2 

<C[||P'«tHia(l|V«|U- +C) + llVdtll^dlVuHL- + ||Vd||^ + C) + HVdtllia + 1] 
<C(\\ P ^u t \\ 2 L 2 + ||Vd t || 2 2)(||Vu|| L oo + \\Vd\\% 2 +C) + C(||Vd t || 2 2 + 1). (2.44) 

Applying the Gronwall's inequality to estimate (I2.44p . we deduce 

sup [ (p\u t \ 2 + \\7d t \ 2 )dx+ [ ||Vu t ||l 2 + ||(Ad-/(d)) t ||| 2 df 
o<t<TJn Jo 
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<C [ (||Vdt||!, + l)dtexp{/ (\\Vu\\ L ~ + \\\7d\\ 2 H2 +C)dt}<C, 
Jo Jo 



(2.45) 



where we have used estimate (|2.15[) and the assumption (|2.ip in the last inequality. This completes 
the proof of Lemma 12.51 n 



The following lemma gives the higher order norm estimates of u, d and p. 
Lemma 2.6 Under the assumption (|2.1[) , it holds that for < T < T* 

sup (\\u\\ H . + \\Vd\\ H 2)<C- (2.46) 

0<t<T 

sup ||Vp|U» + / \\V 2 u\\ 2 L edt<C. (2.47) 

0<t<T Jo 

Proof. From estimates (|2~T4| . ([2~38]) and (f2~43|) . we have 

||u||Ha < C(\\piut\\ L , + \\Vu\\ L 2 + \\Vp\\ L * + HVdtlUO < C. (2.48) 
By using estimates (|2.3ip and (|2.33|) . we have 

l|Vd|| ff 2 <C(||V 3 d|| L 2 + ||Vd|| L 2) 

<C(||Vd t || L 2 + ||V M || L2 + ||M|V 2 d||| L 2 + HVdllia) 

<C(\\Vd t \\ L i + \\Wu\\ L 2 + ||Vn||| 2 + ||Vd|| £a +C)<C, (2.49) 

where we have used the estimates (|2.10|) . (|2 . 14f) and (|2.38|) in the last inequality. Combining the 
estimates (12.481) and (|2.49|) above gives the estimate (|2.46l) . 

Applying the operator V to the mass conservation equation (jTTTJ) , , then multiplying the result- 
ing equation by 6|Vp| Vp and integrating over Q give 

4l|Vp||| 8 =- 6 / \Vp\ 6 Vudx - [ V(|Vp| 6 ) • udx - 6 / |Vp| 6 divwda: - 6 / p\Vp\ 4 VpV div udx 
dt Jo Jo. Jn Jq 

<C||Vu|| L »||Vp||i8 + C||Vdiv«|| z e||Vp||i6 ) 

that is 

^||Vp|Ua <C||V«||zoo||Vp||L»+C||Vdivu||i.. (2.50) 
By using the Gronwall's inequality to the above estimate gives 

||Vp|U« <(\\po\\w^+C [ \\\7divu\\ L edt)exp{C [ ||V«|| L o.di} 



<C(f 
Jo 



T 

2 



V 2 u\\ L 6dt+l). (2.51) 



Applying the standard elliptic regularity result ||V 2 m||l 6 < C||^ u IIl 6 , Holder inequality, Sobolev 
embedding, the estimates (|2.10p and (|2.46[) . we have 

||V 2 U || L6 <C(\\pu t \\ L e + \\pu- Vu\\ L e + \\VP\\ L e + ||(Vd) T (Ad - f(d))\\ L e) 
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<Cr(||Vttt|| £a + ||«|U-||V«|| £ . + HVplU- + ||Vd||L-||Ad|| za + ||Vd|| L a||/(d)|| L ») 

<C(\\Vu t \\ L 2 + \\u\\ 2 H2 + \\Vp\\ L e + \\d\\ 2 H3 + \\d\\ H 2) 

<C(\\Vu t \\ L 2 + \\V P \\ L e + l). (2.52) 
Inserting the estimate (|2.52p into (|2.5ip yields 

||Vp|| z . <C [ {\\Vu t \\ L 2 + \\Vp\\ L e + l)df 
Jo 

<C ( {\\Vu t \\ 2 L 2 + \\Vp\\ L , + l)df < C [ ||Vp|| £ .df + C, 
Jo Jo 

where we have used the estimate (|2.38p . then applying the Gronwall's inequality gives 

sup ||Vp|| L6 < C. (2.53) 

0<t<T 

From (|2"32|) and ([2T53]) . we have 

/ \\S7 2 u\\ 2 Le dt < C( / ||Vw t ||| 3 dt + sup \\Vp\\le +C) < C. 

JO JO 0<t<T 



(2.54) 



It is easy to known that the estimate (|2.48[) follows (|2.53[) and (|2.54l) immediately. This completes 
the proof of Lemma |2~61 □ 



We now give the proof of Theorem 11.21 

Proof of Theorem 11.21 From the existence result of Theorem ll.il we know that ||u(f)||fl-2, 
||p(t)||vpi,e, ||d(t)||jy3 and ||p^itt(i)||£2 are all continuous on the time interval [0,T*). From the 
above Lemmas [2~Tti2T6l we see that for all T € (0,T*), 

(\\u\\ H 2, \\p\\ W i,e, \\d\\ H3 , \\p?u t \\ L 2)(T) < C. (2.55) 

Furthermore, there hold 

p^u t + p^u ■ Vu G L°°([0, T*]; L 2 ), (2.56) 

and for all T G (0,T*), 

(pAu-Adiv(Vd0 Vd-i(|Vd| 2 +F(d) j) -VP)(T) = {pu t + pu ■ Vu)(T) = ^/pg, (2.57) 

where g(T) = (p^ut + p^u ■ Vu)(-,T) G L 2 . Therefore, from (|2.56j) and (|2.57j) . we can take 
(p, u, d)\ t= T with any T G (0,T*) as the initial data and apply Theorem 11.11 to extend the local 
strong solution to a time interval [T, T + 6] for a uniform 6 > which only depends on the bounds 
obtained in these lemmas, so that the solution can be extended to the time interval [0,T* + S). 
This contradicts with the maximality of T*. Hence, the assumption (|2.ip cannot be true. This 
completes the proof of Theorem 11.21 □ 
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